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Abstract 

A study of D + K® and D°K + final states is performed in a sample of l.Ofb -1 of pp 
collision data collected at a centre-of-mass energy of \fs = 7TeV with the LHCb 
detector. We confirm the existence of the L>* 1 (2700) + and D*j(2860) + excited states 
and measure their masses and widths to be 

m(L>* 1 (2700) + ) = 2709.2 ± 1.9(stat) ± 4.5(syst) MeV/c 2 , 

r(Z)* 1 (2700) + ) = 115.8 ± 7.3(stat) ± 12. l(syst) MeV/c 2 , 

m(D* J (2860) + ) = 2866.1 ± l.O(stat) ± 6.3(syst) MeV/c 2 , 

r(D*j(2860) + ) = 69.9 ± 3.2(stat) ± 6.6(syst) MeV/c 2 . 
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1 Introduction 



The spectrum of the known cs states is at present described as two S-wave states (Df, 
D* + ) with spin-parity assignment J p = 0~ , 1~ and four P-wave states (-D* (2317) + , 
L> sl (2460) + , £> sl (2536) + , L>: 2 (2573) + ) with J p = 0+, 1+, 1+ 2+ |1], of which the latter 
two have also been observed in semileptonic -B-decays in LHCb pf. This picture is still 
controversial since the _D* (2317) + and -D s i(2460) + states, discovered in 2003 [3]j6], were 
predicted to have much higher masses [7[jlT]. Between 2006 and 2009, three new D s j 
mesons were observed at the B factories in DK and D*K decay mode^j]and in three-body 



^hadron decays: the D* 1 (2700)+ [^14], the £>* 7 (2860) + pp j] and the L> sJ (3040) + [Mj 
excited states. From the angular analyses in Refs. [l3||l4| , J p = 1~ is favoured for the 
D* 1 (2700) + state, a possible J p = 3~ assignment is discussed for the Z)*j(2860) + , and an 
unnatural parity is suggested for the D s j(3040) + state since it was found to decay only 
to the D*K final state. 

The measured properties of the _D* 1 (2700) + state are in agreement with theoretical 
expectations (7- 10 , but further confirmation is still needed. Similarly, the existence of the 
_D*j(2860) + resonance is unclear. In the latest analysis by the BaBar collaboration [14] , 
the observation of the _D*j(2860) + decaying to the D*K final state rules out the J p = + 



assignment of Refs. 15 , 16 , reinforcing the 3 assignment proposed in Refs. 17 18 . How 



ever, the measured branching fraction ratio £>(.D*j(2860) + -> D*K)/B(D* sJ (2860) + 
DK) = 1.1 ± 0.2 is larger than the predicted values of 0.39 [l7| and 0.6 [18]. In Ref. [19], 
it is argued that the observed decays of the _D*j (2860) + to DK and D*K do not rule out 
a J p = 0+ assignment, but the pattern of decays can be explained by the overlap of two 
radially excited states of J p = + and J p = 2 + , with similar mass and width. 

Given the controversial status of these high mass D s j states, none of them is currently 
reported in the summary table of the Particle Data Group [l] . Experimental contributions 
are needed in order to disentangle the puzzle around the D* sJ (2860) + and to complete the 
picture of the cs spectrum. 

Using 1.0 fb -1 of data recorded by the LHCb detector during 2011 we perform an 
analysis of the D + K® and D°K + final states 

in order to confirm the existence of the D* 1 (2700) + and _D*j(2860) + states and to 
measure their masses and widths. 



2 Detector description 

The LHCb detector [20 is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < r\ < 5, designed for the study of particles containing b or c quarks. The detec- 
tor includes a high precision tracking system consisting of a silicon-strip vertex detector 
surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4Tm, and three stations of silicon- 

1 DK refers to D+K° and D°K+, while D*K refers to D* + K° and D*°K+ final states, where the 
inclusion of charge conjugate final states is implicit everywhere. 
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strip detectors and straw drift-tubes placed downstream. The combined tracking system 
has momentum resolution Ap/p that varies from 0.4% at 5GeV/c to 0.6% at lOOGeV/c, 
and impact parameteij^] resolution of 20 [im for tracks with high transverse momentum 
(pr) with respect to the beam direction. Charged hadrons are identified using two ring- 
imaging Cherenkov detectors. Photon, electron and hadron candidates are identified by 
a calorimeter system consisting of scintillating-pad and pre-shower detectors, an electro- 
magnetic calorimeter and a hadronic calorimeter. Muons are identified by a muon system 
composed of alternating layers of iron and multiwire proportional chambers. The trig- 
ger consists of a hardware stage, based on information from the calorimeter and muon 
systems, followed by a software stage which applies a full event reconstruction. 

Monte Carlo simulated event samples are used to calculate the effects of the detector on 
the mass resolution. The pp collisions are generated using Pythia 6.4 21 with a specific 
LHCb configuration [22]. Decays of hadronic particles are described by EvtGen [23] 
and the interaction of the generated particles with the detector and its response are 



implemented using the Geant4 toolkit 24,25 as described in Ref. |26|. Simulated 



events are reconstructed in the same manner as data. 



3 Event selection 

We reconstruct the D + K® final state using the D + — > K~7t + tt + and K® — > tt + tt~ decay 
modes, and the D°K + final state using the D° — > K~tt + decay mode. Because of their 
long lifetime, K® mesons may decay inside or outside the vertex detector. Those that 
decay within the vertex detector acceptance have a mass resolution about half as large as 
those that decay outside of its acceptance, as observed in Fig. [T] 

Tracks are required to have good track fit quality, momentum p > 3 GeV/ c and trans- 
verse momentum px > 250 MeV/ c. Tracks pointing to a pp collision vertex (primary 
vertex) are rejected by means of an impact parameter requirement in the reconstruction 
of the D + , D° and K® candidates. The tracks used to reconstruct the mesons decaying 
inside the vertex detector are required to have a distance of closest approach among them 
smaller than 0.5 mm. 

To improve the signal to background ratio for the reconstructed D + , D° and K® 
meson candidates, we require the cosine of the angle between the momentum of the 
meson candidate and the direction defined by the positions of the primary and the meson 
decay vertex, to be larger than 0.9999 for K® and 0.99999 for charmed mesons. This 
requirement ensures that the meson candidates are produced in the primary pp interaction, 
and reduces the contribution from particles originating from 6-hadron decays. The D + 
and K®, and similarly D° and K + candidates, are fitted to a common vertex requiring 
X 2 /ndf < 8, where ndf is the number of degrees of freedom. The purity of the charmed 
meson candidates is enhanced by requiring the decay products to be identified by the 
ring-imaging Cherenkov detectors, using the difference in the log-likelihood between the 
kaon and pion hypotheses A In Ck-k- We require AlnCxn > 2(0) for kaon tracks and 

2 The perpendicular distance between the track path and the position of a pp collision. 
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Figure 1: Invariant mass distribution (points) for (a) D + , (b) D°, decaying (c) inside and 
(d) outside the vertex detector. We show the total probability density function (solid curve), the 
signal component as a sum of Gaussian distributions (dotted curve) and a decreasing exponential 
distribution to describe the background component (dashed curve). The region within the 
vertical lines corresponds to ±3 standard deviations in mass resolution from the measured peak. 



AIiiCktt < 10(6) for pion tracks from D + (D°) decays. The overlap region in the particle 
identification definition of a kaon and a pion is small and not a problem given the reduced 
number of multiple candidates per event. Figure [T] shows the invariant mass spectra for 
the D + , D° and K® meson candidates after the described selection is applied. The signal 
regions for D + , D° and K® candidates correspond to ±3 standard deviations in mass 
resolution from the peak values. 

At 7TeV, charged track multiplicities from pp interactions are very high, extending 
beyond 100 tracks per event, leading to large combinatorial background. We define 9 as 
the angle between the momentum direction of the kaon in the DK rest frame and the 
momentum direction of the DK system in the laboratory frame. This variable is symmet- 
rically distributed around zero for resonant states, but more than 90% of combinatorial 
background events are in the negative cos# region. We therefore require cos# > to 
strongly reduce combinatorial background, for both D + K® and D°K + final states. A 
further reduction of this type of background is achieved by performing an optimization 
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Figure 2: Invariant mass distributions for (a) D+K$ and (b) D°K+. 



of the signal significance of the cleanest D s j peak in the DK samples, the _D* 2 (2573) + 
state. In the 2.5 — 2.6GeV/c 2 mass region of the DK spectra, we compute the max- 
imum of the signal significance Ng / v Ns + Nb, where Ns and Nb are the number of 
signal and background events, as a function of different requirements on discriminating 
variables. This study motivates the following choices. For the D + K® final state we re- 
quire p T (D + Kg) > 4.5GeV/c for K® candidates decaying inside the vertex detector, and 
Pt(K®) > 1.5GeV/c for K® candidates decaying outside the vertex detector. For the 
D°K + final state we require pt(K + ) > 1.5GeV/c and Ptssn k (K + ) > 0.45, trained using 
inclusive fully simulated Monte Carlo samples and calculated from a neural network using 
as input particle identification log-likelihoods, momenta, tracking related variables and 



sub-detector acceptance requirements combined with Bayesian statistical methods 27 . 



4 Analysis of the DK invariant mass spectra 

The resulting D + K® and D°K + invariant mass distributions are shown in Fig. [2J where 
we have reconstructed about 0.36 x 10 6 D + K® and 3.15 x 10 6 D°K + candidates with 
a multiplicity of 1.005 and 1.010 candidates per event. The D + K® and D°K + mass 
spectra show very similar features. The sharp peak near the threshold is due to the 
feed-down from D sl (253Q) + D* + K°, D*°K+ decays, with D*+ D+vr , £>+7 and 
D*° D°7T°, D°j, where the neutral pion or photon have not been reconstructed. Since 
the L> s i(2536) + state has J p = 1 + , the decay to DK systems is forbidden by angular 
momentum and parity conservation. The observed feed-down is well isolated and the 
overlap with high mass structures is negligible. A prominent peak is observed around 
2.57GeV/c 2 , corresponding to the spin-2 _D* 2 (2573) + resonance. We also observe two 
broad structures near 2.71GeV/c 2 and 2.86GeV/c 2 in both mass spectra, which previous 
measurements [14] have associated with the spin-1 D* 1 (2700) + state and the D*j(2860) + 
state. 

We perform a binned (5MeV/c 2 bin size) simultaneous extended maximum likelihood 
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fit to the two DK mass spectra in the 2.44 — 3.46GeV/c 2 range, where the lower bound 
excludes the _D s i(2536) + feed-down events. Hereafter we will refer to this as the reference 
fit. 

The D s j signal components are described by relativistic Breit-Wigner lineshapes in- 
cluding the Blatt-Weisskopf form factors which limit the maximum angular momentum 



in a strong decay via the introduction of an effective radial meson potential 28 . Mass 
resolution effects are neglected in the reference fit, since the expected widths for the 
Z)* 1 (2700) + and _D*j(2860) + states are between one and two orders of magnitude larger 
than the detector mass resolution, but these effects are included as a source of systematic 
uncertainty. The background distribution is largely dominated by randomly associated 
DK pairs created during the hadronization processes, and is described using a linear 
combination of Chebyshev polynomials of the first kind, of order from one to six. These 
polynomials are flexible and capable of describing possible background fluctuations from 
non-resonant events. The analytical function to describe the background component was 
trained on a fully combinatorial wrong-sign sample of D°K~ events, reconstructed and 
selected in the same way as the D°K + final state candidates. Additionally, we generate 
a sample of signal events where the D s j components of the probability density function 
are taken from the combined DK and D*K measurement performed by the BaBar ex- 



periment 14 . From the combination of the wrong-sign and signal simulated samples we 
study possible fit instabilities and correlations of the width of the D* 1 (2700) + state as a 
function of the lower fit bound. 

The signal model was chosen from a set of fits to the DK mass spectra, where we 
include and remove the expected D s j states from the fit function, with their masses and 
widths fixed to the previous BaBar measurement. The reference signal model, which shows 
the best x 2 /ndf, includes the spin-2 L>* 2 (2573) + , spin-1 D* 1 (2700)+ and D* sJ (28Q0) + states. 
Regarding the _D s j(2860) + state, we use a spin-0 hypothesis since at present no conclusive 
J p assignment has been made for this state. With the current data sample we are not 
able to identify the presence of additional states in the 2.86 GeV/c 2 region, as proposed in 



Ref. 19 . In order to reduce correlations between the background function and the width 
of the broad resonances and to improve fit stability, we fix the less contributing and most 
correlated parameters, the order three, five, and six Chebyshev polynomial coefficients 
for the two DK invariant mass spectra. These parameters are taken from a preliminary 
fit, where the signal model is fixed to values obtained using an approximate background 
shape, similar to that used in the BaBar analysis 14 and described in Section [5] 

The reference fit includes a total of twenty-six parameters, fourteen to describe the 
background components (six fixed as mentioned above) and twelve for the description of 
the signal contributions. The six parameters for the masses and widths of all the D s j 
structures are constrained to be the same in the D + K® and D°K + samples. The reference 
fit results for the D* 1 (2700) + and _D*j(2860) + parameters and total number of events are 
reported in Table [T] and Table [2j respectively. The projections of the fitted function super- 
imposed to the data and the residuals after subtracting the fitted background distribution, 
are shown in Fig. |3j 

The fit quality is acceptable with a total xV n df of 464/422=1.1. We account for 
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Table 1: Parameters for D* 1 (2700)+ and D* J (2860) + states, evaluated with binned fits to the 
samples. Masses and widths are given in units of MeV/c 2 . Uncertainties are statistical only. 





£>* x (2700)+ 


D;j(2860)+ 


Fit sample 


xViiclf 


m 


r 


m 


r 


Reference fit to D + K® and D°K + 


464/422 


2 709 ± 2 


115 ±7 


2 866 ± 1 


70 ±3 


D+Kl only fit 


207/214 


2 710 ±4 


100 ± 14 


2 867 ±3 


73 ±7 


D°K + only fit 


241/214 


2 709 ± 2 


117 ±8 


2 866 ± 1 


67 ±4 



Table 2: Total number of events for D* 1 (2700) + and D*j(2860) + , evaluated with the reference 
fit. Uncertainties are statistical only. 



Decay mode 


L>: 1 (2700)+ 


£^(2860) + 


D+K° s 


6 724 ± 596 


4 825 ± 347 


D°K + 


45 315 ±2 186 


31 603 ±1257 



imperfections in the magnetic field map and alignment of the tracking system. These cor- 
rections are computed using a sample of D° — > K~tc + decays, using the momentum scale 
calibration method explained in Ref. 



29 . The corrections were found to be compatible 



with zero and therefore neglected. 



5 Cross-checks and systematic uncertainties 

The fit is validated using a large set of simulated experiments. No biases are observed and 
the resolution reported by the fit to data is found to be in agreement with the resolution 
from the analysis of the generated experiments. As a cross-check, we perform a set of 
fits to different data subsamples. We perform independent fits to the D + K® and D°K + 
samples (Table [I]) and to the D + K® sample splitting the contributions from the K® meson 
decaying inside and outside the vertex detector. We repeat the reference fit on different 
DK samples recorded with positive and negative magnet polarity, and also in a data 
sample of candidates required to pass dedicated D + and D° triggers. In all cases, we 
found the fit results to be compatible with the reference fit. 

Systematic uncertainties are summarized in Table [3] There are calculated as the 
difference between the results of alternative fits and the reference fit, unless otherwise 
stated. 

A systematic uncertainty is associated to the signal model. Given the unknown J p 
assignment for the _D*j(2860) + excited state, we repeat the reference fit assuming spin-1, 
spin-2 and spin-3 hypotheses for this resonance. A second systematic contribution to the 
signal description comes from the fact that the Blatt-Weisskopf form factors introduce a 
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Figure 3: Invariant mass distributions (points) for (a) D + K® and (b) D°K + . We show the 
total simultaneous probability density function (solid line), the D* 2 (2573) + (fine dotted line), 
D* 1 (2700)+ (dot-dot-dot dashed line), L>*j(2860) + (dot dashed line) and background contri- 
bution (dashed line). Invariant mass distributions after combinatorial background subtraction 
are shown for (c) D + Kg and (d) D°K + , where the vertical scales are truncated to show the 
L>* 1 (2700)+ and D*j(2860) + signals more clearly. 



penetration radius that we fixed in the reference fit to 1.5 GeV -1 . The contribution to the 
systematic uncertainty is estimated by varying this value within the 1—3 GeV -1 range. 
In both cases, we take the largest variation as systematic uncertainty. The quadratic 
combination of these two effects represents the largest systematic contribution to the 
L>*j(2860) + parameters. 

The background component is highly correlated with the yield and width of the 
broad structures, particularly for the D* 1 (2700) + state. Four uncorrelated effects are 
studied. We use an empirical function to describe the background component in the 
D + K® decay mode. This function, similar to that used in the BaBar analysis 14 , is 
composed of a threshold function multiplied by a decreasing exponential of the form 
(m — m t h) P exp {— c\m — C2fn 2 }, where m t h = m(D + ) + m(Kg). On the D°K + sample, 
this function does not reproduce correctly the background shape. Instead we generate 
a set of samples, using the reference probability density function, but randomly varying 
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Table 3: Systematic uncertainties for the L>* 1 (2700) + and -D*j(2860) + parameters. Mass and 
width uncertainties, 5m and ST, are given in units of MeV/c 2 . The total uncertainties are 
calculated as the quadratic sums of all contributions. 





J D* 1 (2700) + 


23^(2860)+ 


Source 


Sm 


SF 


Sm 


SF 


Signal model 


2.2 


3.0 


5.5 


3.4 


Background model 


2.1 


10.2 


3.8 


4.2 


High mass state 


0.0 


0.3 


0.0 


0.2 


Selection criteria 


2.1 


3.5 


1.0 


2.7 


Mass resolution 


2.1 


3.6 


2.8 


2.4 


Feed-down reflections 


1.2 


2.9 


0.1 


1.4 


Bin size 


0.2 


0.9 


0.0 


0.2 


Total 


4.5 


12.1 


6.3 


6.6 



the background parameters. The average difference between the generated and fitted 
values for the _D* 1 (2700) + and _D*j(2860) + masses and widths is taken as the system- 
atic uncertainty. We repeat the reference fit changing the lower bound of the fit range by 
±10MeV/c 2 and the upper bound by —50 MeV/c 2 . This has the largest effect on the width 
of the _D* 1 (2700) + state since the broad width is sensitive to modifications in the amount 
of background near the threshold and in the long high-mass tail. Finally we evaluate a 
systematic uncertainty given by the effect of fixing some of the background parameters 
in the reference fit. We perform a set of fits accounting for all possible up and down 
variations (independently and simultaneously) of these parameters. The variations are of 
10% for D + K® background parameters and of 5% in the case of the D°K + decay mode. 
According to a fit \ 2 study, alternative fits with larger variations of the fixed parameters 
do not describe the data correctly and therefore not used to compute systematic uncer- 
tainties. We adopt as systematic uncertainty the root-mean-square variation of all the 
fits for the given parameter. As expected, this effect contributes mainly to the widths of 
the resonances since these parameters correlate strongly with the background shape. The 
total background model systematic uncertainty is the quadratic combination of the four 
effects discussed. 

Evidence for an additional broad state around 3 GeV/ c 2 has been shown previously in 



D*K decay modes [14]. Therefore, in addition to the D* 1 (2700) + and D* J (2860)+ high 
mass states, we allow for another signal component in the fit. No statistically significant 
structure is found. 

The uncertainty introduced by the selection criteria is computed by repeating the fit in 
a sample with the following selection: pi(D + K®) > 4.75GeV/c and p^(Kg) > 1.7GeV/c 
for D + Kg combinations with the Kg meson decaying inside and outside the vertex de- 
tector, respectively, while for the D°K + sample we apply p^(K + ) > 1.8GeV/c and 
Pn~n k {K + ) > 0.5. These selection criteria are established by optimizing the signal signif- 
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icance of the _D* 2 (2573) + in the 2.5 — 2.6 GeV/c 2 range, as done previously, but this time 
downscaling the number of signal events by one order of magnitude 0.1Ns/y/0.1Ns + Nb, 
trying to mimic the signal to background ratio observed for the _D* 1 (2700) + and -D*j(2860) 
states. 

Mass resolution effects are neglected in the reference fit since the measured widths 
are much larger than the mass resolution obtained from Monte Carlo simulated 
data: 4.3 (3.3) MeV/c 2 at 2.71 GeV/c 2 and 5.2 (4.0) MeV/c 2 at 2.86 GeV/c 2 mass for the 
D + Kg(D°K + ) decay mode. This effect is accounted for by a convolution of the relativis- 
tic Breit-Wigner lineshapes with a single Gaussian function without offset whose width 
is fixed to the mass resolution estimated using fully simulated events. Here, the largest 
contribution arises from the _D* 2 (2573) + state, since a narrower width for this state causes 
a deviation in the masses and widths of the resonances under study. 

The observed D* 1 (2700) + and D* J (2860) + states can also decay into D*K final states 
(depending on the -D*j (2860) + spin-parity) and this should be reflected as feed-down 
components to the DK samples, arising from D* + — > D + ir°, D +r y and D*° — > D°ir , D°^y 
decays, where the neutral pion and photon are not reconstructed. In this case, we expect 
the feed-down structures to be shifted by about — 142MeV/c 2 from the measured mass 
and with similar width but with a small spread from resolution effects. Ignoring resolution 
effects, we evaluate a systematic uncertainty due to the presence of possible feed-down by 
including the two additional components to describe the D* 1 (2700) + — > D* + K®, D*°K + 
and _D*j(2860) + — > D* + K®, D*°K + processes, with fixed masses and widths to avoid 
large correlations. The uncertainty due to this effect is about a factor two smaller than 
the statistical precision on the masses and widths. 

Finally, to investigate the effect of binning the data samples, we repeat the fit using 
bins with size of 1 MeV/c 2 . This effect is observed to be negligible. 

The total systematic uncertainty is calculated as the quadratic sum of all the men- 
tioned contributions. The systematic uncertainties on the L>* 1 (2700) + and Z)*j(2860) + 
parameters dominate the overall measurement uncertainties. 



6 Conclusions 

Using l.Ofb -1 of data recorded by the LHCb experiment during 2011 in pp collisions at a 
centre-of-mass energy of y/s = 7TeV, we perform a study of the D + K® and D°K + final 
states. We observe for the first time the production of D* 1 (2700) + and _D*j(2860) states 
in hadronic interactions and measure their parameters to be 

m(D* 1 (2700) + ) = 2709.2 ± 1.9(stat) ± 4.5(syst) MeV/c 2 , 

r(D; i (2700) + ) = 115.8 ±7.3(stat)± 12. l(syst) MeV/c 2 , 

m(D* J (2860) + ) = 2866.1 ± l.O(stat) ± 6.3(syst) MeV/c 2 , 

r(Z)* J (2860) + ) = 69.9 ± 3.2(stat) ± 6.6(syst) MeV/c 2 . 



All results are compatible with previous results from the B factories 13p4 . The statistical 



uncertainties for all parameters are improved by an overall factor of two with respect to 



9 



the BaBar measurements in the same decay modes, and it is of the same order as for the 
combined DK and D*K BaBar measurement. The precision of the measured quantities 
is dominated by systematic effects. We do not observe any statistically significant D s j 
resonance in the mass region above 3 GeV/c 2 . 

To shed light on the puzzle around the spin-parity of the -D*j (2860) + state and to 
confirm the spin-parity assignment of the D* 1 (2700) + , an angular analysis of D*K samples 
would be needed. 
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